Tunable adsorption on carbon nanotubes 
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We investigated the adsorption of a single atom, hydrogen and aluminum, on single wall carbon 
nanotubes from first-principles. The adsorption is exothermic, and the associated binding energy 
varies inversely as the radius of the zigzag tube. We found that the adsorption of a single atom and 
related properties can be modified continuously and reversibly by the external radial deformation. 
The binding energy on the high curvature site of the deformed tube increases with increasing radial 
deformation. The effects of curvature and radial deformation depend on the chirality of the tube. 



PACS numbers; 73.22.-f, 68.43.Bc, 68.43.Fg, 68.43.-h 



Novel mechanical, electrical and chemical properties of 
carbon nanotubes ||l|-g|l have been explored actively with 
a motivation of finding a new technological application. 
A single wall carbon nanotube (SWNT) is usually de- 
scribed by a rolled graphene, where the hexagonal 2D 
lattice is mapped on a cylinder of radius R with various 
helicities characterized by a set of two integers (n, m). A 
SWNT can display either metallic or insulating electronic 
structure depending on the hclicity and radius [Ej. 

Recent studies m-m showed that the electronic prop- 
erties of SWNTs can be modified by radial deformation. 
The energy gap of an insulating SWNT can decrease and 
eventually vanish at an insulator-metal transition with 
increasing applied radial strain. The density of states 
at the Fermi energy, 'D{Ef), of a metallized SWNT in- 
creases with the increasing radial strain. More interest- 
ingly, the radial deformation necessary to induce metal- 
licity was found to be in the elastic range. Therefore, 
all strain induced changes in the electronic and also in 
mechanical properties are reversible. 

Most noticeably, the radial strain disturbs the unifor- 
mity of charge distribution. This, in turn, may impose 
changes in the chemical reactivity and hence in the inter- 
action of tube surface with foreign atoms and molecules. 
It is therefore anticipated that not only band gap but also 
chemical reactions taking place on the surface of a SWNT 
can be engineered through radial deformation. In this pa- 
per, we explore this feature by using the predictive power 
of the density functional theory and demonstrate that 
indeed adsorption of foreign atoms on carbon nanotubes 
and associated properties can be modified continuously 
and reversibly. Furthermore, we showed that there is a 
simple scaling of the adsorption energies with the radius 
of the SWNT. We believe that the tunable adsorption 
can have important implications for metal coverage and 
selective adsorption of foreign atoms and molecules on 
the carbon nanotubes, and can lead to a wide variety of 
technological applications, ranging from hydrogen stor- 
age to new materials Q. 

We have investigated the electronic energy structure 



and charge density of bare and single atom adsorbed 
SWNTs with and without radial strain by performing 
first- principle density functional calculations Q. We ex- 
pressed the wave functions by linear combinations of 
plane waves up to an energy cutoff of 500 eV within the 
tetragonal supercell geometry (separation between the 
tubes is taken to be 7.5 A and c is equal to the lattice 
constant of SWNT along its axis). With this energy cut- 
off and using ultra soft pseudopotentials for carbon |l(]| ] 
the total energies of nanotube-adatom systems converge 
within 0.5 meV/atom. The Monkhorst-Pack ||ll[ special 
k-point scheme is used with 0.02 A~^ k-point spacing 
resulting 5 k-points along the tube axis. Results have 
been obtained within the generalized gradient approxi- 
mation |jl3] (GGA) for fully relaxed geometries including 
all carbon and adsorbate positions and the lattice con- 
stant of the tube along the z-axis. 
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FIG. 1. Binding energies Ei, of single hydrogen and alu- 
minum atom adsorbed on the zigzag SWNTs versus the ra- 
dius of the tubes. The solid line is the fit to Eb = Eo -\- C/R 
(see text). 

To investigate the effect of the radial deformation we 
first consider adsorption of H and a simple metal, Al, on 



undeformed SWNTs. In Fig. ||, we present the binding 
energies of H and Al adsorbed on the (n, 0) zigzag nan- 
otubes calculated for n =7,8,9,10,12. H is adsorbed at 
the top site, i.e. directly above the C atoms of the tube, 
Al favors the hollow site i.e. above the center of the 
hexagon as in the graphite surface. The binding energy 
Eb, 

Eb = Et [SWNT] + Ea-Et [SWNT + A] (1) 

is calculated in terms of the total energy of SWNT, 
Et[SWNT], total energy of SWNT with an adsorbed 
atom A, Et[SWNT + A], and the energy of the single, 
free atom, Ea. Here the bare SWNT and SWNT with an 
adatom A are free of external stress, and all the atomic 
coordinates are fully relaxed. Moreover, since E}, is calcu- 
lated by using the same supercell, the spurious adatom- 
adatom interaction along the tube axis is substracted. 
The positive value of Eb indicates that the adsorption is 
exothermic and hence stable. It is found that the binding 
energy of an adatom decreases with increasing radius (or 
decreasing curvature) of the tube, and eventually satu- 
rates at a value corresponding to that on graphene plane. 
Hence the variation of binding energies of H and Al with 
the radius of the zigzag tube fits to the curve given by 
the expression. 



EbAR) = Eo.A + % 



(2) 



where Eq.a is the binding energy of the adatom A (H or 
Al) on the graphene plane. We calculated Eq^h =1.49 eV, 
and Eo^Ai = —0.02 eV. Interestingly, the fitting parame- 
ter Ca is found to be independent of the adatom (H and 
Al) and is equal to '^ 3.14 eVA. Currently, we are extend- 
ing our calculations to see if this relation holds for other 
adsorbates as well. The binding energies calculated for 
(n, 0) SWNTs with n < 8 deviate from the above simple 
scaling perhaps due to the fact that the singlet 7r*-band, 
which is normally in the conduction band, falls into the 
band gap as a result of increased a* — n* mixing at high 
curvature 0,|l3l- Note that, while the band gap shows 
significant change with n, (for example, upon going from 
(8,0) to (9,0) Eg changes from 0.65 eV to 0.09 eV) the 
binding energies vary smoothly with R~^. Increasing Eb 
with decreasing R (or with decreasing n) shows that for 
small R the character of the surface deviates from that 
of the graphene. This finding also suggests that by cre- 
ating regions of different curvature on a single SWNT 
by radial deformation one can attain different values of 
binding energies. 

The radial deformation that we consider in this study 
is generated by applying uniaxial compressive stress on 
a narrow strip on the surface of the SWNT m . In prac- 
tice such a deformation can be realized by pressing the 
tube between two rigid flat surfaces. The radius is de- 
creased in the y-direction, while it is elongated along the 



x-direction. As a result, the circular cross section is dis- 
torted to an elliptical one with major and minor axis, 
a and b respectively. The elliptical radial deformation 
can be described by the magnitude of the applied strain 
and it is defined as eyy = {R — b)/R where R is the ra- 
dius of the undeformed nanotube (with zero strain) . For 
different values of Cyy we carried out full structural op- 
timization under the constraint that the minor axis is 
kept fixed at a preset value by freezing the carbon atoms 
at both ends of the minor axis m. Then, under this 
constraint, the coordinates of the adatom and remain- 
ing carbon atoms and the lattice constant of the tube, c 
are optimized until there remain only forces opposite to 
the applied strain on the fixed atoms, but all other force 
components are less than 0.01 eV/A. 




FIG. 2. (a) Variation of the binding energy Eb of a single 
hydrogen atom adsorbed on a (8, 0) zigzag SWNT as a func- 
tion of the elliptic radial deformation, eyy. The upper curve 
corresponds to H adsorbed on the high curvature region near 
the end of the major axis a {sharp site). The lower curve is 
for the adsorption on the low curvature region at the end of 
the minor axis {flat site) (b) Same as (a) for a single adsorbed 
Al atom. Insets: Ball and stick models and isosurface plot of 
difference charge densities, Ap for each case. (See text). 

Figure eK shows the variation of the binding energy 
Eb of a single hydrogen atom adsorbed on the (8, 0) sur- 



face with the apphed radial strain for two cases: (i) The 
binding energy of a single H atom adsorbed on the high 
curvature side of the surface {i.e. specified as the sharp- 
site near one of the ends of the major axis, a; = a, y = 0) 
traces the upper curve in Fig. 2a. (ii) The lower curve 
is for the adsorption on the low curvature side near one 
of the ends of the minor axis {i.e. x = 0, y = b specified 
as the flat-site). Here, the binding energy is defined as 
in Eq(l), except that Et[SWNT] and Et[SWNT + A] 
are calculated for radially deformed SWNTs. The bind- 
ing energy of H adsorbed on the high curvature (sharp) 
site is increased by 0.85 eV for eyy—0.3. On the other 
hand, Ef, for the adsorption on the low curvature (flat) 
site first decreases with increasing Cyy, and then satu- 
rates at an energy 0.25 eV less than that corresponds to 
eyy =0.0. The difference of binding energies of the sharp 
and flat sites, AEi,, is substantial and is equal to ^ 1.1 
eV. This value is 44 % of the binding energy of H on the 
undeformed SWNT. As a result of H adsorption, the sp^ 
character of the bonding of the tube has changed locally 
and became more like to sp'^. The lengths of the C-C 
bonds at the close proximity of H have increased slightly. 

The binding energy of Al shown in Fig. 0b exhibits 
a behavior similar to that of H, despite H and Al favor 
different sites on the (8,0) tube: Ef, at the sharp site of 
the deformed SWNT increases with increasing eyy. For 
example, Et, increases by ~ 0.80 eV for eyy — 0.3 which 
is 80% of the binding energy on the undeformed tube. 
For Al absorbed on the flat site, Ef, first decreases with 
increasing eyy, then gradually increases. Adsorption of 
Al induces local changes in the atomic and electronic 
structure. For example, the surface of the tube where 
Al is adsorbed expands. 

The variation of Ef, with the radial deformation is con- 
sistent with the results illustrated in Fig. |l|. In gen- 
eral, the higher is the curvature under deformation, the 
higher the binding energy. The effect of the elastic de- 
formation is further investigated by analyzing the charge 
density and electronic structure. The MuUiken analysis 
estimates that ~0.37 electrons is transfered from H to 
carbon [Q. The C-H bond is directional and is cova- 
lent with a partial ionic component. The charge differ- 
ence, Ap = p[SWNT + H]- p[SWNT] - p[H], which 
is calculated in terms of the undeformed (or deformed) 
SWNT with H (adsorbed at different sites), undeformed 
(deformed) clean SWNT, and single H, are shown by in- 
set in Fig. 0a. Ap indicates no significant change with 
radial deformation. On the other hand, the charge trans- 
fer upon the adsorption of Al is different from H. Since 
Al is adsorbed on the center of a hexagon, the bond be- 
tween Al and SWNT is distributed to nearest C atoms 
with charge accumulation between Al and those C atoms. 
The bond charge slightly increases at the sharp site as 
shown by inset in Fig. 0b. The charge transfer from Al 
atom is estimated to be 0.71 electrons. 
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FIG. 3. Electronic band structure, total and partial den- 
sity of states (with dotted lines) of Al atom adsorbed on a 
SWNT. The zero of energy is taken at the Fermi level shown 
by dash-dotted lines. The band due to Al in the band gap is 
shown by dashed line. Al adsorbed: (a) on the undeformed 
(8, 0) tube, (b) on the undeformed (6, 6) tube, (c) on the flat 
site of the (8, 0) tube, (d) on the sharp site of the (8, 0) tube. 
eyy = 0.31 for (c) and (d). 

Explanation of this remarkable and significant change 
of the binding energy with radial deformation is sought 
in the electronic energy structure and the total and par- 
tial density of states. The adsorption of H gives rise to 
a new state which falls in the gap at F and coincides 
with Ep. This state partially overlaps with the conduc- 
tion band when H is adsorbed at the flat site, whereas 
it occurs near the valence band edge at the sharp site. 
A similar situation occurs with the adsorption of Al as 
illustrated in Fig. ||. The band gap is wide open, and a 
resonant state ~ 2.5 eV below Ep and a p-derived ad- 
sorption state which sets Ep falls at the center of the gap 
for the adsorption of Al on the undistorted SWNT [^ . 
In the case of adsorption on the flat site the latter state 
overlaps the conduction band (Fig. ^), but it dips into 
the valence band at the sharp site (Fig. ||d). Accordingly, 
the sharp-site adsorption involves the valence band states 
of SWNT, while the conduction band states are involved 



in the flat site adsorption. 




FIG. 4. The shifts of band edges at F-point with radial 
deformation on a (8, 0) tube. H adsorbed: (a) on the flat site, 
(b) on the sharp site; Al adsorbed: (c) on the flat site, (d) on 
the sharp site. The valence band edge, the conduction band 
edge, the adsorption state and the Fermi energy are shown by 
diamonds, squares, circles and dots, respectively. 

The above argument is in accord with the binding en- 
ergies calculated for the adsorption of Al on the unde- 
formed and deformed (6, 6) armchair SWNT. Note that 
the (6, 6) as well as all (n, n) SWNTs are metallic with fi- 
nite 'D{Ep) owing to the tt~ and tt*— states crossing and 
setting Ep. Our calculations show that 'D{Ep) of the 
(6, 6) decreases slightly with increasing radial deforma- 
tion, but it remains essentially finite. We found Ei, =0.91 
eV for the adsorption of Al on the (6, 6) SWNT. Upon 
the deformation of the tube with Cyy =0.22, Ei, corre- 
sponding to sharp-, and flat-site adsorption has changed 
to 0.95 and 0.85 eV, respectively. Therefore, the differ- 
ence between the binding energies of the sharp and flat 
site adsorption is only 0.1 eV, and hence is negligible as 
compared to the situation above discussed for the (8, 0) 
SWNT. In Fig. Bp the state associated with the adsor- 
tion of Al on the undeformed (6, 6) tube occurs above the 
Fermi level and donates its electron to the metallic tt— 
and TT*— bands. It appears that Ef, is almost pinned by 
the above mechanism. 

Radial deformation induced changes in the energies 
of valance band edge (VBE) and conduction band edge 
(CBE), as well as the adsorbtion state (AE) at F-point 
illustrated in Fig. H corraborate the above arguments. In 
the case of adsorption on the flat site, VBE is almost 
constant for both H and Al, while CBE and AE both 
shift downwards. As a result of this, more conduction 
band states contribute to the bond energy. By contrast, 
for adsorption on the sharp site, while AE behaves very 
similar to the previous case, CBE is now approximately 



constant, but VBE shifts upwards. In this situation, the 
contribution to the bond energy from conduction band 
states is decreased. 

In conclusion, we showed that the chemical reactivity 
of a zigzag SWNT can be modified reversibly and vari- 
ably by the radial deformation. The effect of deforma- 
tion is significantly different for the zigzag and armchair 
SWNTs. It is remarkable that Al, which is not bound 
to the graphite surface, can be adsorbed at a high cur- 
vature site of a zigzag SWNT under radial deformation 
with a binding energy of 1.8 eV. This novel property 
may have important implications for various chemical 
and electronic applications, such as selective absorption 
and desorption of molecules and atoms, fragmentation 
and chemical sensors, magnetic tubes etc. 
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